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ABSTRACT 

The variation in solar heating due to the resonance rotation of 
Mercury produces periodic elastic deformations on the surface of 
the planet. The thermal stress and strain fields under Mercury's 
surface are calculated after certain simplifications. It is shown 
that deformations penetrate to a greater depth than the variation of 
solar heating, and that the thermal strain on the surface of the 
planet pulsates with an amplitude of 4 x 10"^ and a period of 
176 days. 



THERMAL ELASTIC DEFORMATIONS 
OF THE PLANET MERCURY 


INTRODUCTION 

In this paper, we investigate the thermal elastic deformations of the planet 
Mercury by means of equations of equilibrium aiid compatibility in which the 
variation of solar heating due to resonance rotation has been involved. In view 
of the localized nature of the maximum thermal stress and strain induced on the 
surface of Mercury the half-space model is adopted for the analysis. 


TEMPERATURE FIELD 

We consider a half-space in the xyz coordinate system with a plane boundary. 
Let a temperature wave T with cyclical frequency cj and amplitude be propa- 
gated on the plane x = 0 with velocity v in the direction of the y axis. Then the 
temperature wave T is described by 


T = Tq cos CO 



( 1 ) 


Because of the orientation of Mercury relative to the sun during each peri- 
helion passage (Liu, 1971) the variation of solar heating on the surface of the 
planet can be taken into account by the introduction of a factor cos(&;/2v)y. The 
temperature distribution on the surface of Mercury is then given by the expression 


^(o.y.t) = Tq cos (6^t - 2ey) COS ey 


( 2 ) 


where co = 4.2 x lO"’ sec~^ , e = co/2v = 2 x 10”^ cm~^ . 

The problem of determining the temperature field T(x,y,t) is to find the 
bounded solution of the heat-conduction equation 
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for 


0 < X < 00 

-00 < y < 00 

-00 < t 

In equation (3) p is density, c the specific heat and k the thermal conduc- 
tivity. If we assume that Fourier's law holds for the propagation of the temper- 
ature wave, the temperature field is 


^(x.y.t) = T^e-^^cos (cut - Sx - 2ey)cos ey 


(4) 


where 8 = i {oj p c)/{2k) . 

THERMAL STRESS 

On Mercury's surface a progressing temperature wave T(o,y,t) is given in 
a way to account for the variation in solar insolation; it generates the tempera- 
ture field (4) along the equator of the planet. From this temperature field, we 
can find the thermal stress. The components of thermal stress are related to 
a stress function <5 by (Landau and Lifshitz, 1959) 
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The stress components must satisfy the equilibrium conditions 
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and the condition of compatibility 
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where E is Young’s modulus, v Poisson's ratio, and a the coefficient of linear 
thermal expansion. It is noted that the external forces in equation (6) are omitted 
because the variation of body forces is small and simple calculations show that 
their influence is negligible for the case of the planet Mercury. 

Let = 4> -4’’ Then equation (7) is equivalent to a system of two differen- 
tial equations 
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Since the solutions of equation (9) can be included in solutions to the bihar- 
monic equation (8), only a particular integral is needed for the solution of equa- 
tion (8). In addition, the boundary conditions are 
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= 0 and cr I = 0 
0 xy x = o 


( 10 ) 


because there is no vertical and shear stress on the surface x = 0. Thus the 
thermal stress problem is reduced to solving equations (8) and (9) with conditions 
in equation (10). The stress will then be given by equation (5) for $ = 


3 



The periodicity of the temperature variation at the boundary leads to a 
periodic thermal stress under the surface. K we use the local character of the 
solar heating during perihelion passage, the stress function can be expressed 
in the form 


<?!) = kj fj (x) exp (-i ey) 

( 11 ) 

4^ = fjCx) exp (-iey) 

By substituting equation (11) into equations (8) and (9), we obtain ordinary 
differential equations. The unknown coefficients can be determined from the 
boundary conditions in equation (10) by using equation (5). The results of the 
components of thermal stress are 
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These components of stress satisfy the equilibrium conditions as can be 
verified by substituting equations (12), (13) and (14) in equation (6). 
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INCLINATION AND DISPLACEMENT 

The thermal strain and displacement (U,V) are related to thermal 

stress , cr^y) by means of the generalized Hooke’s law 
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— [p + 1 + (y - 1) ex] e”^’' sincJt 
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sin^c^t - 8x +^^j. siney (20) 

LATITUDINAL STRESS AND STRAIN 

The thermal stress and strain are accompamed by a stress component 
given by 
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and by a strain component given by 
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Therefore, the latitudinal stress component a and strain component e are 
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NUMERICAL RESULTS 

The maximum value of the thermal elastic inclination, d, in equation (19) 
occurs at<^t = tt/ 4 for x = 0. This implies that the surface inclination reaches 
its maximum value when Mercury is at a position which is 7 days after peri- 
helion passage. Numerical calculations were carried out for E = 10^' dyne cm" 2 , 
To = 3 X 1Q2 deg, k = 0.5 x 10"2 cal cm"^ sec"^ deg"^ p = 5 gcm"^, c = 0.2 calg"^ 
deg V = 0.3 and a = 10"® deg"*; The main results are: 6 | = 2.6" x 10"^, 

le I =4x10'^ and |(T„„i = 4 x 10® dyne cm" It should be pointed out 

that calculations of equation (19) show that the thermal elastic inclination, 6, 
penetrates to a greater depth than the temperature variation in solar heating 
because it is not only produced by the absolute value of the temperature, but 
also by the temperature gradient. 


CONCLUSION 

The variation of solar heating on Mercury's surface due to resonance rota- 
tion generates periodic thermal elastic inclination and displacement. From the 
foregoing analysis, we may draw the following concluding remarks. 

(1) The thermal stress on the surface of Mercury is about 4 x 10® dyne 
per cm2 

(2) Thermal elastic inclination penetrates to a greater depth than the 
temperature variation of solar heating. The maximum value of sur- 
face inclination occurs when Mercury is in a position which is 7 days 
after perihelion. 

(3) Equations (19) and (20) describe the thermal elastic response to the 
variation of solar heating. The latitudinal strain on the surface of Mer- 
cury's thermal bulges pulsates with an amplitude of 4 x 10*^ in a period 
of 176 days. 
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